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Fluid-enhanced surface diffusion controls 
intraparticle phase transformations 

YiyangLi© 1 - 2 ' 13 , HungruChen 3 , KipilLim 14 , HaitaoD. Deng 1 , JongwooLim 12 , DimitriosFraggedakis© 5 , 
Peter M.Attia® 1 , SangChul Lee 1 , Norman Jin 1 , JozeMoskon 6 , ZixuanGuan 7 , William E.Gent® 8 , 
JihyunHong© 1 - 4 , Young-SangYu 9 , MiranGaberscek 610 , M. Saiful Islam 3 *, Martin Z.Bazant 1,5 ' 11 ' 12 * 
and WilliamC.Chueh 12 * 

Phase transformations driven by compositional change require mass flux across a phase boundary. In some anisotropic solids, 
however, the phase boundary moves along a non-conductive crystallographic direction. One such material is Li x FeP0 4 , an elec¬ 
trode for lithium-ion batteries. With poor bulk ionic transport along the direction of phase separation, it is unclear how lithium 
migrates during phase transformations. Here, we show that lithium migrates along the solid/liquid interface without leaving 
the particle, whereby charge carriers do not cross the double layer. X-ray diffraction and microscopy experiments as well as 
ab initio molecular dynamics simulations show that organic solvent and water molecules promote this surface ion diffusion, 
effectively rendering Li x FeP0 4 a three-dimensional lithium-ion conductor. Phase-field simulations capture the effects of sur¬ 
face diffusion on phase transformation. Lowering surface diffusivity is crucial towards supressing phase separation. This work 
establishes fluid-enhanced surface diffusion as a key dial for tuning phase transformation in anisotropic solids. 


S olid-state phase transformations triggered by the migration of 
atoms and ions enable many important technologies, including 
batteries 1 , hydrogen storage 2 , electrocatalysts 3 and memristors 4 . 
In anisotropic materials, the mass flux required for phase transfor¬ 
mations may lie along a non-conductive crystallographic direction. 
Because bulk mass transport is slow, phase transformations in these 
materials require alternative mass transport mechanisms, which are 
not well understood. 

Li x FeP0 4 (0<X<1) is a model phase-transforming anisotropic 
material 5-8 with a large miscibility gap for particles larger than -50 nm 
(ref. 9 ) and exhibits one-dimensional lithium migration 10,11 . At low (de) 
lithiation rates, Li x FeP0 4 intercalates via a particle-by-particle pathway, 
first explained by the domino cascade model 12 and more recently by its 
non-monotonic Li chemical potential profile 13 (that is, mosaic instabil¬ 
ity 14 ). At elevated (de)lithiation rates, density functional theory 15 and 
reaction-limited phase-field modelling 16 suggest that phase separation 
is suppressed and replaced with a solid solution pathway. This meta¬ 
stable solid solution was observed using operando diffraction and 
microscopy 17-21 , and the solid solution phase separates into Li-rich 
and Li-poor phases under equilibrium at room temperature 15 (Fig. la). 
Owing to anisotropic strain, the stable phase boundaries are perpen¬ 
dicular to the [100] and [001] directions 14,21-23 , indicating that lithium 
traverses along these directions during phase separation. Yet lithium- 
ion diffusion is fast in the [010] direction and essentially non-existent 
in the other directions 10,11,24 . These two central features of Li x FeP0 4 
present an important but unresolved question: if the bulk crystal struc¬ 
ture does not transport lithium in the [100] and [001] directions, then 
how does lithium migrate during phase transformations? 


There are several possible avenues for lithium transport along 
nominally non-conductive directions (Fig. lb). Anti-site defects 
enable ion hopping between the fast conduction channels 25-28 , but 
this bulk resistance R hu]k is very large for particles with low anti¬ 
site defect concentrations 11,25 . Surfaces may enable higher trans¬ 
port rates, but density functional theory calculations show that 
the lithium surface diffusion resistance R sm{D at the FeP0 4 /vacuum 
interface is even higher than P bulk (ref. 29 ). Since in-plane bulk and 
surface diffusion are believed to be slow, the prevailing model is 
that lithium de-intercalates from the particle, migrates in the liquid 
electrolyte and re-intercalates 16,18,19,30 . Here, the interfacial reaction 
resistance P rxn is assumed to be rate-limiting because of fast lithium 
ion transport in the electrolyte 31 . Many-particle phase-field simu¬ 
lations show that this pathway cannot result in intraparticle phase 
separation because the electrode prefers to phase separate between 
particles 32,33 . Such simulations, however, contradict in situ observa¬ 
tions of intraparticle phase separation within individual particles 20 , 
and highlight an insufficient understanding of the lithium migra¬ 
tion pathway and the phase separation mechanisms in Li x FeP0 4 . 

Phase transformation in ionically isolated particles 

To elucidate the phase transformation mechanisms, we first mea¬ 
sured the rate of phase separation in ionically isolated, solid solu¬ 
tion Li 05 FePO 4 microplatelet particles 21,23 (Supplementary Figs. 
1-3). Both carbon-coated and uncoated particles were studied. 
These particles are -3 pm in the major axis, which is nearly par¬ 
allel to [001], ~1.5pm in the minor axis, nearly parallel to [100], 
and -170nm thick in the [010] thickness direction (Supplementary 
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Fig. 11 Crystallographic directions, phase separation and in-plane lithium 
migration in a Li 05 FePO 4 platelet particle, a, Solid solution Li x FeP0 4 
separates into Li-rich and Li-poor phases. The phase boundaries lie 
along the ab and be planes, and perpendicular to the [100] and [001] 
directions. Therefore, lithium must migrate in the [100] and [001] 
directions. The [001] direction is the major axis along the plane, whereas 
the [100] direction is the minor axis, b, Cross-section schematic view of 
the crystallographic directions. Three possible in-plane migration paths 
are possible: bulk diffusion, surface diffusion and electrolyte diffusion in 
conjunction with interfacial (de)lithiation reactions. The dashed green 
lines indicate that electrolyte diffusion enables lithium transport between 
particles. The relative resistances of these three paths (R bu , k , R surfD , and R rxn ) 
dictate the lithium migration path taken during the phase transformation. 


Figs. 1 and 4). While larger than technologically relevant battery 
particles, this particle morphology has been extensively used in 
fundamental studies 21-23 , and is thin enough to minimize the effect 
of blocked conduction channels 25 . We estimate the anti-site defect 
concentration to be between 0.1-0.5% based on Rietveld refinement 
(Supplementary Table 1). 

Figure 2a shows the evolution of the X-ray diffraction (XRD) pat¬ 
terns as the electro chemically prepared, carbon-coated Li 05 FePO 4 
electrode relaxes in an inert Ar glovebox (<lppm 0 2 and H 2 0) at 
room temperature (see Methods for sample preparation details). 
We observe the expected (020) peaks for the Li-rich triphylite (LFP) 
and Li-poor heterosite (FP) at 29.8° and 30.9°, respectively. Because 
this electrode was electrochemically lithiated at 2 C, where C/N is 
the rate to (de)lithiate in N h, it exhibits significant solid solution 
behaviour, with the strongest peak at 30.1° likely to be a metastable 
phase 17,34 . The (211) peaks are weak due to preferential orientation 
(Supplementary Fig. 4). As the electrode relaxes over hundreds 
of hours, the solid solution peaks diminish while the LFP and FP 
peaks grow, confirming that this electrode slowly phase separates. 
To quantify the solid solution fraction (SSF), we fit Gaussians to 
the diffraction patterns (Supplementary Fig. 5 and Supplementary 
Section 3) 19 . The quantified SSF decreases from 57% to 33% after 
500 h in Ar. In Li 05 FePO 4 particles without carbon coating (Fig. 2b), 
the SSF also decreases slowly, from 39% to 26% in 100 h. Both elec¬ 
trodes retain significant solid solution character over long times, 
reflecting a high R bulk . 


Next, we examined the effect of solvent exposure on the rate of 
phase transformation in uncoated particles, electrochemically lithi- 
tated at 1C to Li 0 5 FeP0 4 . When the electrode is immersed in a salt- 
free ethylene carbonate/dimethyl carbonate (EC/DMC) solvent, 
the SSF decreases from 41% to below 30% in just 8h, more than 
one order of magnitude enhancement in the phase separation rate. 
Ambient atmosphere similarly increases the phase separation rate. 
Solvent, protons 35 , oxygen and water molecules do not intercalate 
into the bulk of Li x FeP0 4 , in part due to large differences in size from 
the ionic radii of lithium. XRD and transmission electron micros¬ 
copy detected no secondary phase formation on the (010) facet 
from exposure to ambient water vapour (Supplementary Fig. 7). 
This is consistent with the stability of Li x FeP0 4 with regard to 
both oxygen and hydrogen evolution, and its use in aqueous Li-ion 
batteries 36 . Although we cannot fully rule out chemical reactions 
between water vapour and the Li x FeP0 4 surface, which occur for 
extended exposure to liquid water 37 , the dominant effect of water 
vapour here is to facilitate phase separation of the underlying solid 
solution Li x FeP0 4 . Because the phase separation rate is propor¬ 
tional to the lithium migration rate 38 , we propose that solvent and 
water molecules interact with the lithium on the surface of the par¬ 
ticles and activate surface diffusion. 

We repeated this experiment for carbon-coated particles, lithi¬ 
ated at 2C (Fig. 2c). Exposing Li 05 FePO 4 to water vapour or to 
ambient air increases the phase separation rate by at least two orders 
of magnitude: 8h of exposure to llTorr of water vapour in Ar 
induces significantly more phase separation than 500 h of exposure 
to pure Ar (Fig. 2c). The salt-free EC/DMC solvent does not accel¬ 
erate phase separation in the carbon-coated particles, in contrast 
to the uncoated particles. The role of the carbon coating in light of 
our experimental results is discussed in Supplementary Section 3 
and warrants further work. We later discuss phase separation in the 
presence of the electrolyte. 

We also directly visualized the distribution of lithium within 
individual carbon-coated particles using scanning transmission 
X-ray microscopy (STXM, Fig. 2d). The lithium composition X is 
computed from the Fe oxidation state maps by fitting a linear com¬ 
bination 21,32 of LiFe 2+ P0 4 and Fe 3+ P0 4 references (Supplementary 
Fig. 8a) to the acquired absorption spectra at every pixel. The vis¬ 
ible Li-rich, solid solution, and Li-poor domains are significantly 
larger than the pixel size of 50 nm. According to the STXM results, 
the particles in the electrode stored in inert Ar are mostly solid 
solution, whereas the ones stored under ambient air are mostly 
phase-separated. The SSF, defined as the fraction of pixels where 
0.15 <X< 0.85, is 84 ±12% in Ar and only 32 ±10% in air. STXM 
and XRD here provide complementary information: microscopy 
confirms that the environment directly alters the lithium compo¬ 
sition within particles, and that the observed diffraction patterns 
do not arise only from elastic strain. Likewise, the significant SSF 
computed from the (020) reflection (Fig. 2a) confirms that mixed 
valence (yellow) pixels in STXM do not result from a phase bound¬ 
ary parallel to the (010) facet, or from phase-separated domains too 
small to be captured with our ~50nm pixel size. Supplementary 
Figure 9 captures relaxation within a single particle and confirms 
that phase separation results from spinodal decomposition rather 
than nucleation and growth. 

Atomistic origins of enhanced surface diffusion 

Having established that water and solvent molecules enhance the 
rate of phase separation in Li 05 FePO 4 by at least two orders of 
magnitude, we investigated the atomic-scale mechanism of lith¬ 
ium diffusion at the solid/fluid interface. We performed ab initio 
molecular dynamics (MD) calculations on LiFePO 4 (010)/water, 
LiFePO 4 (010)/EC and bare LiFePO 4 (010)/vacuum interfaces, build¬ 
ing on our previous surface simulation work 39 . Figure 3a illustrates 
the simulation snapshots indicating the migration of a surface Li + 
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Fig. 2 | Tracking solid solution Li 05 FePO 4 electrodes during phase separation, a, X-ray diffraction of carbon-coated particles as they phase-separate in 
inert Ar. The (020) peaks for the Li-rich triphylite (LFP) at 29.8° and Li-poor heterosite (FP) at 30.9° grow over time, whereas the aggregate intensity of the 
intermediate peaks, quantified by the solid solution fraction (SSF), diminish. t = 0 corresponds to the XRD patterns ~10 min after the current is interrupted, 
b, The phase separation rate for bare, uncoated particles is heavily increased in solvent and ambient environments. The error of 0.01 equals the difference 
in the quantified SSF on consecutive scans of the same sample, c, Phase separation for carbon-coated particles is also strongly affected by the environment. 
The filled squares represent in situ diffraction of an electrode relaxing in the battery electrolyte in Fig. 4a. Each marker colour and shape combination in 
b and c represents a different electrode; lines are to guide the eye. d, Direct STXM imaging of the lithium distribution via Fe oxidation state maps confirms 
that air significantly enhances the rate of phase separation compared with Ar. The [001] direction is the projected major axis of each platelet; the [100] 
direction is the minor axis; the [010] direction is the thickness of the platelet and perpendicular to the plane of the page. The thickness in the [010] 
direction is labelled next to each particle, n equals the number of imaged particles, and the STXM SSF confidence interval equals two standard errors. 


at the LiFeP0 4 (010)/EC interface. The migrating lithium is located 
initially at a surface site at the top of a [010] channel. As the simula¬ 
tion proceeds, the Li + coordinates to two EC molecules and a lat¬ 
tice oxygen. The EC molecules lift the Li + from its original position, 
and then move the Li + to a surface site at the top of a neighbour¬ 
ing [010] channel about 4.6 A away. Interestingly, the EC molecules 
are involved in the migration process rather than acting as simple 
adsorbates, and do not intercalate into the bulk of the material. 
Other migration pathways are also possible, and the virtue of our 
MD approach is that such information on molecular and lithium- 
ion dynamics would not be extracted from purely static simula¬ 
tions based on nudged elastic band methodology. In any case, our 
ab initio results clearly indicate that the EC molecules facilitate the 
surface migration of Li + . While not directly simulated, it has been 
shown that electrons migrate via polaron hopping in the bulk 26 . 

MD calculations were also performed for a LiFePO 4 (010)/ 
water interface. In a similar fashion to the LiFePO 4 (010)/EC inter¬ 
face, water molecules were found to coordinate to Li + , assisting its 
migration across the surface, illustrated as simulation snapshots in 
Fig. 3b. In contrast, migration of Li + is not found on the 
LiFePO 4 (010)/vacuum interface during the simulation timescale, 
which is consistent with the high lithium migration energy (~0.9 eV) 
that we calculate on the bare inert (010) surface and agrees with 
a migration barrier previously reported 29 . Moreover, this result is 
consistent with the extremely slow rate of phase separation when 
the particles are stored in Ar gas (Fig. 2). The surface structure con¬ 
sidered for these simulations may be modified with coatings and 


other phosphate phases 40,41 , which warrants further investigation. In 
general, the ab initio MD simulation results indicate that the fluid 
molecules significantly enhance the surface diffusivity, and that this 
solvent-mediated mechanism is the dominant pathway for in-plane 
lithium-ion migration. 

Phase transformation in ionically connected particles 

Next, we investigated phase transformation in ionically connected, 
carbon-coated particles by relaxing solid solution particles towards 
equilibrium in a liquid battery electrolyte (1M LiC10 4 in EC/DMC). 
Here, all lithium redistribution pathways are available (Fig. lb). 
In situ and ex situ X-ray diffraction confirm that the particles phase- 
separate at a rate much faster in the presence of the liquid electrolyte 
than in inert Ar (Figs. 2c and 4a). Significant phase separation is 
observed within 2 h of relaxation. 

If surface diffusion is faster than the insertion reaction 
C^surfD < ^rxn> Fig- 4b), then lithium is confined to its original parti¬ 
cle. Consequently, each solid solution particle initially separates into 
Li-rich and Li-poor domains 14 , conserving the amount of lithium 
in each particle. On the other hand, if surface diffusion is slower 
than the insertion reaction (R surfD > R rxn ) and lithium redistribution 
occurs via the electrolyte (Fig. 4c), then lithium is also free to move 
to nearby particles under fast electrolyte transport 31 . Here, the par¬ 
ticles directly relax into the globally lowest energy configuration, 
a mosaic of Li-rich and Li-poor particles that minimizes interfa¬ 
cial 15,42 and strain 14,43 energies associated with intraparticle phase 
separation. This mosaic pattern, also known as interparticle phase 
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Fig. 3 | Atomic geometries at the LiFeP0 4 (010)/EC and LiFeP0 4 (010)/H 2 0 interfaces from ab initio MD simulations indicating Li-ion migration at the 
surface, a, LiFeP0 4 (010)/EC interface, b, LiFeP0 4 (010)/H 2 0 interface. Top and bottom panels show views from the side and above, respectively, at three 
simulation snapshots with the migrating Li + circled by a dashed line. In the views from above, only the outermost surface Li + ions are shown for clarity. We 
note that the (010) surface structure is intrinsically Li-deficient relative to the bulk lattice, containing 50% vacant Li sites in the surface layer 39 . Key: Li ions 
(green); 0 ions in LiFeP0 4 and EC (red); 0 in water (blue); H in water and EC (white); C in EC (black); P0 4 (purple); Fe0 6 (brown). 


separation, arises during (dis)charge at low rates 13,44 . While the elec¬ 
trode will always eventually relax into the lowest-energy mosaic 
(Fig. 4b, c), an initial observation of intraparticle phase separation 
would demonstrate that surface diffusion is the dominant pathway 
for lithium redistribution. 

To test this hypothesis, we used STXM to image the lithium dis¬ 
tribution within particles. _We define actively relaxing particles as 
the particles where 0.15 <X <0.85; here, X is the average lithium 
composition within a particle. Figure 4d-f shows the lithium dis¬ 
tribution in particles from Li 05 FePO 4 electrodes lithiated at 2C 
and relaxed under different conditions. In the control experiment 
(Fig. 4d), the electrode is stored in inert Ar, and the particles remain 
as solid solutions (SSF 75 ±8%). The active particle fraction here 
is 81 ± 14% (Fig. 4d). Figure 4e shows a different part of the same 
electrode that is relaxed for ~8h in the electrolyte. The active par¬ 
ticle fraction decreases slightly to 71 ± 16%, but the SSF decreases 
significantly to 41 ±8%. This result demonstrates that lithium is 
mostly confined to its original particle during phase separation. 
After relaxing for 500 h in electrolyte, the active particle fraction 
decreases to 38 ±20% (Fig. 4f). 

The observation that individual particles initially separate into 
Li-rich and Li-poor domains confirms that lithium is mostly con¬ 
fined to its original particle in the first few hours. Phase-field simula¬ 
tions of two relaxing solid solution particles (Supplementary Fig. 10, 
described in detail in the next section) verify that this sequence 
is only observed when surface diffusion is faster than interfacial 
reaction. The fast surface diffusion during relaxation is a result of 
lithium not having to cross the double layer or undergoing charge 
transfer, and is consistent with our results for uncoated Li 0 5 FeP0 4 in 
EC/DMC (Fig. 2b) and our MD simulations (Fig. 3a). 

Finally, we show evidence of phase separation in the cur¬ 
rent-voltage response of the battery cell. We probed the elec¬ 
trochemical resistance of a solid solution electrode as it relaxed 
at open-circuit (see Supplementary Section 5b for details). Our 
results clearly show that the resistance of Li 05 FePO 4 increases 
over time across a variety of conditions (Supplementary Fig. 11), 
doubling over the course of 50 h. In contrast, a solid solution 
Li 07 (Ni 1/3 Co 1/3 Mn 1/3 )O 2 electrode exhibits negligible change in the 
resistance over time (Supplementary Fig. lib). Phase separation 
increases the electrochemical resistance due to concentration of 
current in a smaller number of active particles 32 , decreased active 
area for intercalation area, and/or lower exchange current density 
for the phase-separated lithium compositions 21 . The tens of hours 


required to increase the electrochemical resistance is consistent 
with the long times required for intraparticle and/or interparticle 
phase separation (Fig. 4). 

Phase-field simulations of surface diffusion 

We used continuum modelling to explore the effect of surface diffu¬ 
sion on the phase separation behaviour of Li x FeP0 4 under constant 
current conditions. The existing Allen-Cahn reaction model with 
generalized Butler-Volmer kinetics 30 predicts suppression of phase 
separation by electro-autocatalysis upon lithium insertion 14,16,21,45 . 
Here, we unify the Allen-Cahn and Cahn-Hilliard models by add¬ 
ing an in-plane diffusion term to simulate the surface diffusion 
demonstrated experimentally in the previous sections. The result¬ 
ing generalized reaction-diffusion equation is given by: 


ax 

dt 


v 


r£*(i -x)vn Li 

k B T 


+ r(X,// Li ,// ext ) 


( 1 ) 


where X is the spatially resolved lithium composition, t the time, D eff 
the average in-plane diffusivity, which is proportional to the sur¬ 
face diffusivity D surface , [i u the local chemical potential of lithium in 
Li x FeP0 4 , [i ext = yw Li + + /u e - the external lithium chemical potential 
established by the electrons in the current collector and the lithium 
ions in the electrolyte, k B the Boltzmann constant, and T the tem¬ 
perature. These terms control the local rate of lithium insertion r, 
which is proportional to the current density j, and governed by 
generalized Butler-Volmer kinetics 30 described in Supplementary 
equation (4) and Supplementary Section 6. 

We analysed the linear stability of equation (1) using the for¬ 
malism derived by Bazant 45 , which determines whether a small 
perturbation in the lithium composition from a homogeneous 
solid solution X=0.5 will grow (unstable) or decay (stable) (see 
Supplementary equations (5) and (6), and Supplementary Fig. 12 
for definitions). Linear instability is a necessary but insufficient 
prerequisite for phase separation. The neutral stability threshold j c 
depends on the current density, the surface diffusivity D surface and 
the exchange current density j (Fig. 5a). As D surface -> 0, the neutral 
stability threshold is proportional to j and results from the auto- 
inhibitory form of the composition-dependent j Q (refs 14 > 16 > 21 > 45 ). We 
denote this as regime I. Increasing D surface transitions into regime II, 
where the neutral stability threshold is proportional to D surface . 
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No electrolyte exposure 
SSF: 0.75 ± 0.08 
Active: 0.81 ± 0.14 


X=0 


Li y FePO A 


8 h electrolyte 
SSF: 0.41 ± 0.08 
Active: 0.71 ± 0.16 


X= 1 
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500 h electrolyte 
SSF: 0.23 ± 0.98 
Active: 0.38 ± 0.20 




Solid solution 


Phase separation 
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Fig. 4 | The distribution of lithium in ionically connected particles shows that phase separation occurs primarily via surface diffusion, a, In situ XRD 
patterns show that solid solution particles separate into Li-rich and Li-poor phases during relaxation; SSF quantification is shown in Fig. 2c. b,c, Schematic 
illustrations of intraparticle and interparticle phase separation, d-f, STXM images of the lithium composition distribution under different electrolyte 
exposure times show solid solution under no electrolyte exposure (d), intraparticle phase separation after short, 8-h electrolyte exposure (e) and 
interparticle mosaic under long electrolyte exposure (f). The active particle fraction is the fraction of particles with an average composition X between 0.15 
and 0.85. The confidence interval represents two standard deviation under a binomial distribution for the active particle fraction or two standard error in 
terms of the SSF. The particles have the same crystal orientations as in Fig. 2d, and the thickness in the [010] direction is labelled next to each particle. An 
unpaired t test of the active particle fraction between e and f yields t = 2.6 and a p value of 0.01. 


To determine if the perturbations grow fast enough to induce 
phase separation, we simulated a particle lithiating at constant cur¬ 
rent as a function of the applied current density j partic i e and D surface 
using equation (1). Figure 5b plots the final computed SSF, defined 
as the fraction of the simulated particle where 0.15 <X< 0.85 when 
X = 0.5. In regime I on the left side of the plot, we observe three 
regions: linearly stable solid solution (blue), linearly unstable quasi’ 
solid solution (blue) and phase separation (black), consistent with 
previous work 14,16 . Supplementary Figure 12 illustrates the differ¬ 
ences between the three regions. The threshold current density j 
approximately defines the boundary between the quasi-solid solu¬ 
tion and phase separation behaviour. 

As we increase D surface in regime II, the threshold j becomes 
proportional to T> surface . When j <j , surface diffusion redis¬ 
tributes lithium from a solid solution into Li-rich and Li-poor 
domains by transporting lithium to the moving phase bound¬ 
ary (Fig. 5c). On the other hand, when j >j , there is not enough 
time for surface diffusion to redistribute the lithium, creating 


(quasi) solid solution particles. Here, lithium incorporates 
into the particle near its original desolvation site (Fig. 5d). 
In both cases, the entire (010) surface participates in the reac¬ 
tion. Based on experimental observations of phase separa¬ 
tion at a particle rate above 0.1 C (Supplementary Fig. 13), 
our simulations suggest that D surface is at least 10 _12 cm 2 s _1 , signifi¬ 
cantly higher than the 10 15 cm 2 s _1 calculated for a bare FeP0 4 / 
vacuum interface 29 . This quasi solid solution appears even when 
the shape of j 0 (X) does not yield complete suppression of phase 
separation by electro-autocatalysis (Supplementary Fig. 14). 

Our simulations show that decreasing D surface slows phase separa¬ 
tion driven by electro-autocatalysis to give quasi-solid solutions in 
regime II. When surface diffusion is too slow to redistribute the lithium 
(Fig. 5d), the particle behaviour will be governed by reaction kinetics; 
as a result, electro-autocatalysis with asymmetric j (X) with respect to 
X = 0.5 results in abroken symmetry in the growth rate ofnon-uniform 
domains between lithium insertion and extraction 45 , as observed in 
previous experiments 21 . 
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Fig. 5 | Origin of phase separation and solid solution from single-particle simulations, a, The analytical linear stability condition as functions of the 
exchange current density and surface diffusivity. The calculations were conducted at X =0.5. Instability in regime I is driven by reaction, whereas 
instability in regime II is driven by surface diffusion, b, The simulated solid solution fraction (SSF) atX =0.5. j is defined as the SSF = 0.5 contour and 
separates solid solution from phase separation, j plots the linear stability curve from a. The exchange current density used was 0.005 Am -2 , previously 
measured 21 , c, During phase separation where i partide < j s in regime II, lithium desolvates on the entire surface, but lithium migrates on the surface before 
incorporating into the bulk at the phase boundary, d, During solid solution where j partic | e > j 5 , lithium must be incorporated into the particle near the original 
desolvation and charge transfer site. Simulated results of phase separation and solid solution in 300 nm particles are shown below the schematics, where 
the colours have the same meaning as in Fig. 4. 


Supplementary Section 7 and Supplementary Fig. 16 extend 
this discussion to a many-particle electrode. The kinetic phase 
diagrams’ in Fig. 5 and Supplementary Fig. 16 establish how elec- 
tro-autocatalytic reaction kinetics, Crate and surface diffusivity 
conspire to control the phase transformation pathway for Li x FeP0 4 , 
in both single and many particles. Our simulation suggests that sur¬ 
face diffusion should be minimized to obtain (quasi) solid solutions 
within individual particles (Fig. 5b and Supplementary Fig. 16c), 
and that Crate should be high relative to the exchange current den¬ 
sity to obtain concurrent, non-mosaic intercalation in the porous 
electrode (Supplementary Fig. 16). 

Our demonstration of fluid-enhanced surface diffusion in 
electrochemically active materials introduces new paradigms 
for controlling phase transformations in anisotropic mixed ionic 
and electronic conductors, where transport could be dominated 
by surface diffusion. We show that engineering the solid/fluid 


interface controls the bulk phase transformation behaviour. 
Moreover, because solid solution particles have lower interfacial 
resistances (Supplementary Fig. 11), suppressing surface diffu¬ 
sion and phase separation may yield higher roundtrip energy 
efficiencies. Minimizing surface diffusion without altering the 
bulk solid could be achieved by either modifying the fluid (for 
example, introducing small amounts of electrolyte additives) or 
by coating the surface. Surface diffusion may also be suppressed 
altogether in solid-state batteries without electrolyte solvents. 
We recognize, however, that phase separation is an intrinsic 
thermodynamic property of the material, and that even dramati¬ 
cally decreased surface diffusion will not completely suppress 
phase separation because of the parallel electrolyte and bulk dif¬ 
fusion pathways (Fig. lb). We also note that lowering surface 
diffusivity may simultaneously create undesirable effects such as 
reduced reaction rate. 


NATURE MATERIALS | www.nature.com/naturematerials 









































































NATURE MATERIALS 


ARTICLES 


In a broader context, the ability to suppress phase separation by 
minimizing surface diffusion is crucial towards lowering coherence 
strain, reducing mechanical fracture and creating longer-lasting 
devices 46 . Such insights may be especially relevant for high-volume- 
expansion electrodes such as silicon. Beyond batteries, the control 
of phase transformation through surface diffusion is also applica¬ 
ble to anisotropic geometries such as thin films, layered materials, 
nano wires, and highly porous structures where the bulk diffusion 
rates vary significantly in different directions. For example, mois¬ 
ture-activated surface diffusion tunes memristive switching in 
thin films 4 . Surface diffusion can control phase transformations in 
alloying silicon nanowires 47 , with one possible pathway schemati¬ 
cally shown in Supplementary Fig. 17. Layered materials includ¬ 
ing graphite, transition metal dichalcogenides, and MXenes often 
exhibit anisotropic diffusivities and phase separation on ion inser¬ 
tion 48 ’ 49 ; if bulk diffusivities between or within layers are slow, sur¬ 
face diffusion may provide important contributions towards phase 
transformation kinetics. Our generalized reaction-diffusion model 
can quantitatively simulate this behaviour across diverse materials 
systems. Finally, by activating transport at the solid/liquid interface, 
nominally insulating materials may now be used as solid electrodes, 
electrolytes and electrocatalysts, such as the surface diffusion of 
lithium polysulfide species in the Li-S battery 50 . 

In summary, we have demonstrated that surface diffusion of lith¬ 
ium ions facilitated by solvent molecules makes Li x FeP0 4 a ‘three- 
dimensional’ conductor. This surface diffusion is a missing link that 
controls the phase transformation rate in Li x FeP0 4 and determines 
the current threshold between phase separation and solid solution. 
This current threshold is well within the range of operation for typi¬ 
cal Li-ion batteries, such as the time-averaged discharge of portable 
electronic devices over one day. We show that the balance between 
insertion kinetics, cycling rate and the surface ion diffusivity gov¬ 
erns the phase behaviour of Li x FeP0 4 . By establishing that select 
molecules can be used to tune the surface diffusion rate, our work 
demonstrates a facile method to control phase separation, which 
is especially relevant for anisotropic and nanosized materials. We 
show that surface diffusion must also be considered for both mech¬ 
anistic understanding and improved device performance. More 
broadly, we provide a general framework for phase transformation 
mechanisms in anisotropic materials where the phase boundary 
movement and ion flux directions are orthogonal. 
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Methods 

Synthesis of LiFeP0 4 particles. LiFeP0 4 microplatelets were solvothermally 
synthesized using a mixed water and polyethylene glycol solvent based 
on our previous work 21 . All precursors for solvothermal synthesis were 
purchased from Sigma-Aldrich. 6 ml of 1M H 3 P0 4(aq) was mixed with 24 ml of 
polyethylene glycol 400. Afterwards, 18 ml of 1M LiOH (aq) was added to create the 
creamy-white Li 3 P0 4 precipitate. This mixture was bubbled with dry N 2 for ~16h in 
a Schlenk line to remove dissolved oxygen. 6 mmol of FeS0 4 -7H 2 0 was dried under 
vacuum in a Schlenk line for ~16h. Next, we transferred 12 ml of deoxygenated 
H 2 0 to dissolve the dried FeS0 4 powder. The lime-green FeS0 4 solution was 
transferred without air exposure to the Li 3 P0 4 suspension, and the entire mixture 
was transferred to a 100 ml Teflon-lined autoclave in a glovebag. The autoclave was 
heated to 140 °C for 1 h, then to 210 °C for 17h, and cooled. The majority of the 
particles reported in this manuscript are ~170nm thick and came from the same 
synthesis batch. The particles used in Figs. 2d and 4a, and Supplementary Fig. 13 
came from a different synthesis batch and are -250 nm thick due to batch-to-batch 
variations, and show very similar quantified SSF when relaxed in Ar according to 
Figs. 2d and 4d. 

After the synthesis completed, the white LiFeP0 4 particles were centrifuged 
three times with deionized water and dried. Carbon-coating was conducted by 
mixing the LiFeP0 4 with sucrose at a mass ratio of 5:1 (LiFeP0 4 :sucrose) without 
breaking the primary particles. This sample was heated to 600 °C for 5 h in a tube 
furnace under flowing Ar to yield the carbon-coated LiFeP0 4 . Bare, uncoated 
particles were annealed under the same condition without mixing with sucrose. 

Next, we assembled the electrodes containing 65 wt% LiFeP0 4 particles (coated 
or uncoated), 25 wt% C65 carbon black and 10wt% polyvinylidene fluoride 
binder (MTI). The suspension was cast onto carbon-coated A1 foil using a doctor 
blade with a nominal thickness of 10 pm, and dried under vacuum at 90 °C for 
-12 h. Next, coin cells were assembled in an Ar glove box using a 9-mm-diameter 
LiFeP0 4 electrode, 2 Celgard separators and a 1.2-mm Li metal counter electrode 
(Alfa Aesar). The electrolyte was 1M LiC10 4 dissolved in a 1:1 (wt/wt) mixture of 
ethylene carbonate (EC) and dimethyl carbonate (DMC). To measure the capacity, 
the electrodes were cycled at 15mAg -1 (Supplementary Fig. 3) between 2 V and 
4 V. After holding the electrode voltage at 4 V for 6 h on the completion of constant 
current cycling, the remaining LFP (020/211) peak is -8% of the combined LFP 
and FP (020/211) peaks 4 intensities, suggesting that 8% of the particles are inactive 
or disconnected (Supplementary Fig. 6). This explains most of the reduced capacity 
of our electrodes (150mAhg -1 obtained versus 170 mAh g 1 theoretical). Using 
STXM, the quantified lithium composition of the particle used as the LiFeP0 4 and 
FeP0 4 reference are 0.97 and 0.02, respectively, and are spatially homogeneous 
(Supplementary Fig. 8). 

Preparation of samples for ex situ XRD and STXM. To prepare the Li 0 5 FePO 4 
electrodes with significant initial solid solution fractions, we electrochemically 
inserted half a formula unit of lithium into a delithiated FeP0 4 composite electrode 
at a constant current of 1C (for uncoated particles) or 2 C (for carbon-coated 
particles) at room temperature (Supplementary Fig. 3). C/N describes the current 
to (de)lithiate the electrode in N h. The electrolyte is 1.0 M LiC10 4 dissolved in a 
1:1 mixture of EC and DMC. Next, we extracted the electrodes from the battery 
cells, rinsed them with excess DMC solvent to remove the electrolyte, dried the 
electrodes under vacuum and stored the electrodes in different environments, 
including inert Ar, air, salt-free EC/DMC mixture, and electrolyte. 

X-ray diffraction. Except for the synchrotron powder diffraction used in 
Supplementary Fig. 2, X-ray powder diffraction was conducted using a Cu 
K-a source and a one-dimensional detector in a Bragg-Brentano geometry. 

We used both a Bruker D8 Advance and PANAnalytical X’Pert. Ex situ powder 
diffraction was conducted on the electrodes after they are extracted from 
the coin cells and rinsed with excess DMC. Because the diffractometer is not 
located in an inert glovebox, the electrode was exposed to ambient atmosphere 
for about 10 min for each ex situ measurement before being placed back into the 
appropriate environment. 

In situ X-ray powder diffraction consisted of a battery half-cell housed in a 
plastic pouch. The pouch cell design is based on our previous work 51 , but the 
plastic pouch housing used here contains an -2 pm aluminium lining to further 
reduce water and oxygen diffusion through the pouch. The Li x FeP0 4 positive 


electrode, lithium foil negative electrode, Celgard separator, battery electrolyte, and 
nickel and aluminium current collector tabs were assembled into the pouch. The 
electrochemical cycling profile of this pouch cell is identical to that of coin cells. 

During cycling, the pouch was clamped between two metal plates. After 
lithiating the positive electrode at 2 C for 15 min, we removed the pouch cell from 
the clamp and conducted X-ray diffraction on the electrode sitting in the pouch. 
This removal process took less than 5 min, and the sealed pouch was not opened 
at any point during the measurement. X-ray diffraction was conducted in a Bruker 
D8 in a standard reflection mode geometry using a 6-26 scan between 26 values 
of 28.5° and 32°. 50 measurements taking -12 min each were conducted. The solid 
solution fraction for the in situ measurements in Fig. 4a were quantified identically 
to those of the ex situ measurements in Fig. 2a-c. The Li x FeP0 4 electrode remained 
inside the electrolyte environment of the cell throughout the in situ measurement, 
and was not exposed to air at any point. 

Scanning transmission X-ray microscopy. The extracted electrodes were 
sonicated in dimethyl carbonate and dispersed onto lacey carbon transmission 
electron microscopy (TEM) grids, and placed in the X-ray microscope. The 
sample was stored in a glovebox and sealed in an aluminium-lined plastic 
pouch for transportation. STXM was conducted at the Advanced Light 
Source, beam line 5.3.2.2 at the Fe L 3 edge. Supplementary Section 4 and 
Supplementary Fig. 8 contain more details on the STXM analysis protocol 
including the thickness determination. 

Ab initio calculations. Ab initio calculations were based on density functional 
theory as applied successfully to a range of lithium battery materials 52-54 . Here, 
the VASP code 55 was used with the PBEsol exchange-correlation functional 56 
and an energy cut-off for plane-waves of 400 eV. Interactions between core and 
valence electrons are described by the projector augmented wave (PAW) method. 
The gamma-point is used for the Brillouin zone sampling. Molecular dynamics 
simulations were performed using the NVT ensemble with a Nose-Hoover 
thermostat. Each interface is simulated at 500 K with a time step of 1.5 fs for 15 ps. 
The LiFeP0 4 (010)/H 2 O interface is modelled by a seven-layer LiFeP0 4 slab along 
the b axis and a cross-section of a 1X 2 ac plane, and 8 water molecules above the 
surface. The LiFeP0 4 (010)/EC interface is modelled using the same LiFeP0 4 slab 
configuration and 4 EC molecules above the surface. 

Supplementary methods. The data analysis and the phase-field modelling 
protocols are given in the Supplementary Information. 

Code availability. The code for the phase-field model as well as the ab initio 
calculations will be made available upon reasonable request to the authors. 

Data availability 

All experimental data within the article and its Supplementary Information will be 
made available upon reasonable request to the authors. 
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